Strong diamagnetic response of metamaterials 
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We demonstrate that there is a strong diamagnetic response of metamaterials, consisting of open or closed 
split ring resonators (SRRs). Detailed numerical work shows that for densely packed SRRs the magnetic per- 
meability, fJ,(u>), does not approach unity, as expected for frequencies lower and higher than the resonance 
frequency, ujq. Below u>o, A t ( tJ ) gives values ranging from 0.9 to 0.6 depending of the width of the metallic 
ring, while above Uq, fj,(u) is close to 0.5. Closed rings have /i ~ 0.5 over a wide frequency range indepen- 
dently of the width of the ring. A simple model that uses the inner and outer current loop of the SRRs can easily 
explain theoretically this strong diamagnetic response, which can be used in magnetic levitation. 

PACS numbers: 75.20.-g, 41.20.Jb, 42.70.Qs, 73.20.Mf 
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Several types of regular materials exhibit diamagnetic be- 
havior characterized by an induced magnetic moment oppo- 
site to the external magnetic field so that their magnetic sus- 
ceptibility x 1S negative. Most of the natural diamagnetic ma- 
terials have very low values of x, of the order of 10~ 6 . The 
largest known diamagnetic value of a natural material is that 
of pyrolytic graphite which is equal to —4.5 x 1CP 4 . A strong 
diamagnetic response is very important since it can be used 
in magnetic levitation. There are two cases of strong diamag- 
netic response: One is that of the superconducting state which 
shows \ — — 1 ( or weaker in the mixed state) practically over 
all frequencies. The other is that of magnetic metamaterials to 
be examined in this work. 

Over the last seven years artificial materials have been 
designed and fabricated consisting of an assembly of units 
(sometimes called magnetic "atoms") which exhibit a resonant 
response to electromagnetic field of appropriate polarization 
driving thus the magnetic permeability /i = 1 + x to negative 
values for a narrow frequency range just above the resonance 
frequencyi 2 ^. Given the fact that the magnetic response is 
in general a weak v 2 /c 2 effect 4 -, it is not unreasonable to ex- 
pect that fi would approach unity away from the resonance. 
However, as it can be seen from Fig. 1, this is not the case. 
Unusually large values of |x| appear leading to values of ft, as 
low as 0.5. More specifically, there is a frequency regime be- 
low the resonance (Fig. la) extending over at least two orders 
of magnitude where p: stays clearly lower than one; the larger 
the width w of the ring (see caption of Fig. 1), the smaller 
the value of fj, is and the wider this frequency regime is. On 
the other side of the resonance, there is again a plateau where 
fi has a value of about 0.5 for the chosen values of the pa- 
rameters. This value is independent of the width w, and it 
is practically equal to the value of fi appearing if the ring is 
closed (See Fig. lb showing a constant /j, clearly less than one 
over a very broad frequency range extending over 3^ orders of 
magnitude). We point out that the diamagnetic strength shown 
in Fig. lb or la is comparable to that of the superconducting 
state. 

While the resonance region has been studied exten- 
sively, the other diamagnetic regimes have received little 
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FIG. 1: (Color online) Re fi vs. log 10 (a/A) obtained by simula- 
tions for a periodic assembly of open (a) or closed (b) square metallic 
rings. The corresponding Im fi vs. log 10 (a/A) are shown in panels 
(c) and (d), respectively. The simulations have been done using the 
commercial software packet Comsol MultiPhysics, which is essen- 
tially a finite-element method frequency domain solver. The complex 
current distribution was obtained directly from the simulated field 
distributions and, averaged over the unit cell, gave the magnetization 
(This averaging produces unphysical results when the wavelength A 
becomes smaller then the unit cell size as evidenced by Im n < 0.). 
Using (4) together with the incident field Ho we found the effective 
permeability shown. Unit cell size is ax ax a± = 10 x 10 x 2 /im 3 . 
Size of the ring is£x£ = 8x8 fim 2 ; its cross-section is d x w, 
where the depth d — 400nm and the width w takes the values 1/im 
(red, solid curve), 2/im (green, dashed curve), and 3/im (blue, dot- 
ted curve). The metal is Drude-model gold with experimental values 
for the plasma frequency, u v = 2n x 2184 THz, and collision fre- 
quency, u T — 2iv x 6.5 THz 9 . 



attention^ 7 -. In this letter we focus on these regimes and we 
explain their main unexpected features. We start with the ap- 
proach of Gorkunov et al. & , according to which its inductance 
L, its capacitance C, and its resistance R characterize each 
open ring. For closed rings there is no capacitance; the mu- 
tual inductances L nm among the rings (which are periodically 
placed to form a lattice) are also taken into account. The final 
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FIG. 2: (Color online) Re /i vs. log 10 (a/A) according to Eq. (1) 
for open rings (a) for the values of the width w = 1/xm (red, solid 
curve), 2fim (green, dashed curve), and 3/im (blue, dotted curve) 
and for closed rings (b) of any width. The corresponding Im \i 
vs. log 10 (a/A) is shown in (c) and (d), respectively. All parameters 
are as in Fig. 1. 



result for /j, and for open rings is 
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where B = fi nA 2 Q 2 / L, n is the concentration of rings 
n = l/(a 2 a±), A is the area of each ring, A s» (£ " ~ 
Q 2 /L = (L eff + ^ nA 2 )-\ L cS = L + £ m L nm 
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= 1/VLC, and 7o = R/L. For 
closed rings f2 m = 0. Notice that [i^n A 2 = (fioA/a±)f 1 
where fioA/a± is L e g in the solenoid limit, and f x = A/ a 2 is 
the two dimensional filling factor. Thus in the limit of close- 
packed rings 
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where f 2 is smaller than but close to one. 

In Fig. 2 we plot /i vs. log 10 (a/A) according to Eq. (1). 
We have taken into account that R and L (to a lesser degree) 
depend on ui because of the skin effect 4 . We see that the reso- 
nance region as well as the plateau above is reproduced fairly 
well, although the resonance is stronger and sharper than in 
Fig. la. Similarly, the closed ring result is in good agreement 
with the simulations data of Fig. lb except of the final high 
frequency rise towards one, which we attribute to the fact that 
at this high frequency the wavelength A m in the medium is 
comparable to twice the unit cell size; then the uniform field 
assumption, on which the retrieval of // from the simulation 
and Eq. (1) is based, breaks downifi. 

On the other hand, Eq. (1) fails to reproduce the diamag- 
netic response below the resonance for the open ring case. To 
understand the reason behind this failure, we return to the sim- 
ulations and we analyze the current distribution as shown in 
Fig. 3. In case (a), below the resonance, the opposite flowing 
currents at the two edges induce a magnetic field which can- 
cels the external magnetic field inside the metallic wire. The 



wider the wire, the more extensive the diamagnetic volume is 
and the lower /i is. At resonance (case (b)), the strong clock- 
wise current at the inner edge induces a strong field which in 
the inner area of the ring cancels or dominates over the exter- 
nal field producing thus the possibility of negative /i. Finally, 
in the plateau above the resonance, the current flows clock- 
wise near the outer edge inducing a field which cancels the 
external field over the area enclosed by the outer edge of the 
ring; this is the reason that the value of fi in this regime is in- 
dependent of the metallic width w (assuming constant £) and 
equal to the closed ring value, since in the latter case the cur- 
rent flow is as in case (c) over the entire diamagnetic regime. 
The asymmetry of the circular current in case (c) is due to the 
superposition of growing linear currents caused by a electric 
cut-wires resonance at higher frequency. 

To model this more complicated behavior of the open ring, 
we introduce a two loop representation of it sharing the single 
capacitance of the gap as in Fig. 4. By defining the effective 
inductances L\=L\ + J2 m ^i,nm. L 2 = L 2 + Y.m L 2, nm , 
L\i = L 2 i = L 12 + J2 m L i2,nm, we find using Kirchhoff's 
circuit equations the following matrix relation connecting the 
currents ij and I 2 in each loop to the external magnetic field 
Ho 



Z\\ Z\ 2 

Z 2 i z 22 



II I A l 

iujuoHq I A 



(3) 



where Zu = —icuLi — (iu>C) 1 + Ri, i = 1,2, Zi 2 = Z 2 \ — 
—iujLi 2 — (iujC)^ 1 , and A\ (A 2 ) is the area enclosed by the 
outer (inner) edge of the ring. The permeability /i is given, 
following Ref. 6 by 
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where the average magnetization M is equal to n(AiIi + 
A 2 I 2 ). In what follows we shall employ the solenoid ap- 
proximation (which is good for close packed rings, i.e. small 
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FIG. 3: (Color online) Current distribution in the SRR ring (the 
color/gray level is the azimuthal current density, the arrows indicate 
current direction at the edges) for the three regimes shown in Fig. 1 
when the external magnetic field is normal to the plane of the rings 
and pointing out of it: (a) In the diamagnetic regime below resonance 
the current is confined near the outer and the inner edges flowing in 
opposite directions enclosing the metal filled area of the ring trace; 

(b) at resonance the current is mainly confined near the inner surface; 

(c) above resonance the current flows clockwise and essentially at the 
outer edges of the metallic ring. 
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FIG. 4: Equivalent circuit for an open ring. Loop 1 (2) represents the 
current path roughly one skin depth wide along the outer (inner) edge 
of the ring. The two loops have a mutual inductance L12 = L21. 



2 
1.5 
I 

0.5 




-4 -3 -2 

login ( a / A) 



1 



2r 

1.5 - 

-3 !- 
^0.5 - 

£ 0- 
-0.5 - 
-1 



5 -4 



-3 -2 -1 

logio ( a M) 




,(a/A) 



■3 -2 -1 

login 



FIG. 5: (Color online) (a) Re \i vs. log 10 (a/A) for the open SRR 
rings with w = 1/im (red, solid curve) together with the contribu- 
tion of inner (blue, dotted curve) and the outer (green, dashed curve) 
edges, (b) shows the corresponding Im fi vs. log 10 (a/A) (red, solid 
curve) together with the contribution of inner (blue, dotted curve) 
and the outer (green, dashed curve) edges, (c) and (d) show Re fi 
vs. log 10 (a/A) and Im n vs. log 10 (a/A), respectively, for the three 
cases discussed in Fig. la with ring widths 1/im (red, solid curve), 
2/im (green, dashed curve), and 3/im (blue, dotted curve). These 
results were obtained according to Eqs (3) and (4). 



aj_) in order to avoid the tedious lattice sums involved in the 
definition of L\, L\, and L\ 2 , which just renormalize pa- 
rameters but do not qualitatively alter the physical behavior. 
According to this approximation we have L\ ps /j,qAi/cl±, 
L 2 ~ n A 2 /aj_, and L 12 = L21, £12 ~ fM)A 2 /a±. 

In Figs. 5a,b we plot fi(uj) together with the correspond- 
ing contributions from the outer and inner loops, while in 
Figs. 5c,d we plot n(w) for the three cases shown in Fig. la 
according to Eqs. (3) and (4). We see from Fig. 5a that the 
diamagnetic contribution below resonance is coming from the 
outer loop, which prevails over the contribution of the inner 
loop, which is paramagnetic. At resonance the inner loop 
dominates, while the outer loop give a negative peak there. 
Finally, above resonance the outer loop gives a constant dia- 
magnetic contribution, while the inner loop contribution ap- 
proaches smoothly the value \i — 1. The Im (Fig. 5b) 



confirms the predominance of the inner loop at the resonance 
frequency. The contribution of the inner loop to Im fi(uj) in 
the diamagnetic region below resonance is negative, while the 
total Im n(co) is, of course, positive. This indicates that the 
outer loop feeds the inner loop with energy there. For compar- 
ison, we calculated the total fj,(uj) within the two-loop model 
for open rings and for the three cases shown in Fig. la. The 
result is shown in Fig. 5c and it is in good agreement with the 
simulations presented in Fig. la. We calculated also, within 
the two loop model, fi(u>) for closed rings and we found that 
it is constant and negative and equal to the high frequency 
asymptotic value for open rings; this result is not plotted, since 
it practically coincides with those in Figs, lb and 2b. Thus the 
two loop model reproduces the /1 vs. to dependence, as de- 
termined by Femlab simulation, as well as all the features of 
the current distributions in the various regimes, and provides 
a clear explanation for the complicated diamagnetic behavior 
of the open rings. 

In conclusion, the present work shows that there are five 
frequency regions, in the magnetic response of open rings. 
This behavior can be understood in term of the outer and 
the inner current loop of each ring proposed here: The low 
frequency regime, uo < R/L, exhibits practically no mag- 
netic response, fj, = 1, because of the resistive damping as 
expected for most mat erials. The second frequency regime, 
R/L < lo < l/V LC, exhibits an unusually strong diamag- 
netic behavior (even for metal volume filling ratios as low 
as 10%) which is due to the outer current loop, while the 
inner one makes a paramagnetic contribution (the two cur- 
rents are connected along the gap); the net result depends on 
the difference A\ — A 2 of the are as en closed by these two 
loops. The third regime, to ps 1 /V LC, is the resonance re- 
gion, which is dominated by the inner loop; this presents a de- 
sign challenge, since what creates a strong diamagnetic back- 
ground below the resonance frequency, namely a small area 
inner loop, makes the resonance weak. The fourth regime, 
1 / V LC < lu < ire/ a, shows again a very strong diamagnetic 
response and depends only on the outer loop, being propor- 
tional to Ax. Finally, in the regime ire/ a < us, the assump- 
tion of the wavelength X m in the medium being much larger 
than the lattice constant a breaks down and both the averaging 
procedure employed in the simulations and the description in 
terms of a simple effective electric circuit fail. 
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